Superfluid Density in a Highly Underdoped YBa 2 Cu306+y Superconductor 
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The superfluid density p B (T) = l/X 2 (T) has been measured at 2.64 GHz in highly underdoped 
YBa2Cu306+i;, at 37 dopings with T c between 3 K and 17 K. Within limits set by the transition 
width AT C w 0.4 K, p s (T) shows no evidence of critical fluctuations as T — > T c , with a mean- 
field-like transition and no indication of vortex unbinding. Instead, we propose that p s displays the 
behaviour expected for a quantum phase transition in the (3 + l)-dimensional XY universality class, 
with p s0 oc (p - po), T c oc (p - p c ) 1/2 and p s {T) oc (T c - T) 1 as T -> T c . 



PACS numbers: 74.72.Bk, 74.25.Nf, 74.25.Bt, 74.25. Ha 

Current research on high temperature superconductiv- 
ity focuses on the underdoped cuprates, in a region of the 
phase diagram where d-wave superconductivity gives way 
to antiferromagnetism [TJ. One proposal for this regime 
is that at temperatures up to about 100 K above T c , 
superconductivity persists locally, with long-range phase 
coherence suppressed by fluctuations in the phase of the 
superconducting order parameter [2J [31 |H [H]. 
Early results showing a linear relation between T c and 
the superfluid density p s (T = 0) [H] provided the origi- 
nal motivation for this point of view, suggesting that T c 
is low in underdoped materials because the phase stiff- 
ness is low. Further support for this idea has come from 
measurements showing a finite phase stiffness above T c 
at terahertz frequencies [10] . and from experiments that 
appear to detect the phase-slip voltage of thermally dif- 
fusing vortices in the normal state [11] . If the physics of 
the underdoped cuprates is indeed that of a fluctuating 
d-wave superconductor there should be a regime where 
quantum fluctuations come into play as T c falls to zero 
with decreasing doping. Here we test this idea with a de- 
tailed study of the doping dependence of the superfluid 
density in the vicinity of the critical doping for supercon- 
ductivity. 

High homogeneity of T c is particularly difficult to 
achieve in the underdoped cuprates, where the control 
parameter is chemical doping and the materials are well 
away from plateaus or turning points in T c (y). The 
YBa2Cu306+y system has two advantages in this doping 
range: with careful work, there can be sufficient control 
of doping homogeneity to produce samples with sharp 
superconducting transitions [TH [T3] ; and the process of 
CuO-chain ordering can be harnessed to provide contin- 
uous tunability of the carrier density in a single sample, 
with no change in cation disorder [13]. This is possi- 
ble because the loosely held chain oxygen atoms in these 
high quality samples remain mobile at room temperature 
and gradual ordering into CuO-chain structures slowly 
pulls electrons from the C11O2 planes, smoothly increas- 
ing hole doping over time [TSl EE]- For this experiment, 
single crystals of YBa2Cu306+j/ were grown in barium 



zirconate crucibles and have high purity and low defect 
levels, with cation disorder at the 10 -4 level [T7]- A 
crystal 0.3 mm thick was cut and polished with AI2O3 
abrasive into an ellipsoid of revolution about the sample 
c-axis, 0.35 mm in diameter. The oxygen content of the 
ellipsoid was adjusted to 06. 333 by annealing at 914°C in 
flowing oxygen, followed by a homogcnization anneal in 
a sealed quartz ampoule at 570°C and a quench to 0°C. 
At this point the sample was nonsuperconducting. Af- 
ter allowing chain oxygen order to develop at room tem- 
perature for three weeks, T c was 3 K. The sample was 
then further annealed at room temperature for six weeks 
under hydrostatic pressure of ~ 30 kbar, raising T c to 
17 K. The sample was cooled to — 5°C, removed from the 
pressure cell, and then stored at — 10° C to prevent the 
oxygen order from relaxing. All further manipulation of 
the ellipsoid was carried out in a refrigerated glove box at 
temperatures less than — 5°C. In between measurements 
of surface impedance, periods of controlled in- situ an- 
nealing at room temperature and ambient pressure were 
used to generate a sequence of 36 dopings as T c relaxed 
back to 3 K. We emphasize that no oxygen entered or 
left the sample during these annealing steps. Subsequent 
rcanncaling under hydrostatic pressure, for a further six 
weeks, returned T c to the starting value of 17 K, where 
the sample was remeasured to demonstrate the reversibil- 
ity of the technique. 

Measurements of the afr-plane surface impedance Z s = 
R s + iX s were carried out at 2.64 GHz by cavity pertur- 
bation, using a sapphire hot-finger to position the sample 
at the i/-field antinode of the TE015 mode of a rutile di- 
electric resonator [13]. All data sets reported here were 
taken with the c-axis of the ellipsoid oriented along the 
microwave magnetic field H r f to induce afo-plane screen- 
ing currents. The surface impedance of the sample has 
been obtained from the measured cavity response using 
the cavity perturbation formula A/b(T) — 2iA/ (T) = 
r(i? s + iAX s ), where Afg(T) is the change in band- 
width of the TE015 mode upon inserting the sample into 
the cavity, and A/o(T) is the shift in resonant frequency 
upon warming the sample from base temperature to T. 
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FIG. 1: (color online). a6-plane surface impedance at 
2.64 GHz for the YBa2Cu306+i/ ellipsoid, at two dopings. 
Ra(T) is measured directly in the experiment. Absolute re- 
actance is obtained by offsetting AX S (T) so that R s and X s 
match in the normal state. Inset: a fluctuation peak in <7i(T) 
at one doping. AT C is set to the difference between inflection 
points in ai (T) on opposite sides of the transition. 



r is a scale factor that applies to the data set as a whole, 
and is empirically determined using a Pb-Sn replica sam- 
ple to an accuracy of 2.5%. The absolute surface reac- 
tance is set in the usual way by matching R s and X s in 
the normal state, where we expect the imaginary part 
of the microwave conductivity to be very small. This is 
illustrated in Fig. [TJ which shows the surface impedance 
Z s = R s + iX s of the ellipsoid at two of the dopings. The 
rounded shoulders in Z S (T) are the result of fluctuations. 
This is seen more clearly in the microwave conductivity 
a — <j\ — i(T 2 , which is obtained using the local- limit ex- 
pression a = luipLo/Z 2 . The inset of Fig. [l] shows &i(T) at 
one of the dopings, revealing a narrow fluctuation peak 
at T c . In a homogeneous system, u\{T) is expected to 
have a sharp cusp at T c [19J. Rounding of the fluctuation 
peak, arising from the macroscopic dopant inhomogene- 
ity present in all real samples, is used to define AT C . 

The superfluid density is given by p s = 1/A 2 = top^u^- 
Figure [2] shows p s (T) at 20 of the 37 dopings. The 
most striking feature of the data is the wide range of 
linear temperature dependence, extending from close to 
T c down to T w 4 K. Below 4 K p s (T) crosses over to 
an accurately quadratic temperature dependence. Such 
behaviour is well established in YBa 2 Cu306+i, at higher 
dopings and is consistent with <i-wave superconductivity 
in the presence of a small density of pair-breaking de- 
fects [2U] [2T]. One unusual feature of the data is the 
nature of the thermal transitions, which, within limits 
set by AT C appear mean- field-like. At optimal doping 
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FIG. 2: (color online). afr-plane superfluid density p s (T) = 
1/X 2 (T) shown at 20 of the 37 dopings measured in this study. 
Measurements were made starting in the most ordered state 
(T c m 17 K) followed by controlled oxygen annealing in small 
steps down to T c ~ 3 K. At the end of the experiment the sam- 
ple was reordered and measured again to verify reproducibil- 
ity. Lines mark where the vortex-unbinding transition should 
occur for a 2D superconductor. The dashed line corresponds 
to pf° = h 2 d/4k B e 2 p X 2 = (2/it)T in each CuO a plane. The 
solid line shows p 2I> = (2/ty)T in each CUO2 bilayer. p B (T) in- 
stead passes smoothly through this region. While mean-field- 
like over most of the doping range, p s (T) develops downwards 
curvature near T c at the highest dopings, a possible indication 
of the onset of the 3D-XY critical fluctuations. 



YBa2Cu30g+j, is the most three dimensional cuprate, 
with \ 2 C {T 0)/X 2 ab (T -» 0) sa 50 |23|. Its criti- 
cal behaviour has been firmly established to be in the 
3D-XY universality class [11 Hfl US]. In the doping 
range explored in this paper, YBa2Cu306+t, is highly 
anisotropic, with X 2 C (T -> 0)/X 2 ab (T -> 0) w 10,000 [TJ]. 
In these circumstances, one would anticipate fluctuations 
in adjacent layers to be uncorrelated, and a Kosterlitz- 
Thouless-Berezinsky (KTB) vortex unbinding transition 
[HI [2S] should occur when the 2D phase stiffness in a 
layer of thickness d, pf{T) = h 2 d/4k B e 2 p a X 2 (T), falls 
to (2/tv)T. This defines minimum superfluid densities 
for isolated planes and CuC>2 bilayers, shown in Fig. [2] 
by the dashed and solid lines respectively. p s (T) instead 
passes smoothly through these lines, with no indication 
of vortex unbinding. Surprisingly, this implies that fluc- 
tuations remain correlated over many unit cells in the c 
direction. Recent work on YBa 2 Cu306 +y thin films sup- 
ports this, showing that the KTB transition does occur 
but that the effective thickness for fluctuations is the film 
thickness [27] . At the highest dopings in our experiment, 
p s (T) develops slight downward curvature near T c , pos- 
sibly indicating the emergence of 3D-XY criticality. 
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FIG. 3: (color online). Our p s o{p) data (solid diamonds, 
right-hand scale) and T c (p) data from Ref. [13] (solid circles, 
left-hand scale). p s o = 1/A 2 (T — > 0) is obtained by linear 
extrapolation to T = 0. Hole doping for the superfluid density 
data is determined from our values of T c , using a linear fit and 
extrapolation of the T c (p) data from Ref. [T3] (solid line). To 
the extent that this extrapolation is valid, Pso(p) ~ {p — Pc) 2 
at the onset of superconductivity, before crossing over to a 
linear doping dependence. The dashed line is a linear fit to 
p s o(p) at higher doping. 

Previous studies of superfluid density in the under- 
doped cuprates have focused on the strong correlation 
between T c and p s0 = 1/A 2 (T -> 0) HH]. Our new 
data, taken on a single, high purity crystal, explore the 
highly underdoped region in considerably more detail. 
The strong correlation between T c and p s o remains, but 
PsO is up to an order of magnitude larger than in the 
earlier work on thin films |28j . As in the thin film 
study, p s o falls continuously to zero on underdoping and 
varies approximately quadratically with T c at low dop- 
ing. Recent experiments on the CuC>2 plane doping state 
in YBa2Cu30g+j, have established a mapping between 
T c (y) and p, the hole concentration per planar Cu |13j . 
In Fig. [3j a linear fit to the T c (p) data from Ref. US] has 
been used to determine p and plot p s o (p) . At higher dop- 
ing pso has a linear doping dependence. This behaviour 
appears to be very robust: an extrapolation of the linear 
fit in Fig. [3] passes within 5% of the a6-averaged super- 
fluid density of Ortho-II YBa 2 Cu 3 6+a [22 . To the ex- 
tent that the linear extrapolation of T c (p) holds, p s o(p) 
varies approximately quadratically close to the onset of 
superconductivity. However, as we will discuss below, 
this quadratic behaviour is difficult to understand theo- 
retically. Later on we will present an alternative proposal 
in which it is p s o, not T c , whose linear doping dependence 
extends to the edge of the superconducting phase. 

The very low superfluid density, its continuous vari- 
ation with doping, and the importance of ZD-XY fluc- 
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FIG. 4: (color online). The results of assuming a linear 
mapping between doping and p s o (solid diamonds). T c (open 
squares) varies as (p — Pc) 1 ^ 2 as p — > p c , with the solid line 
a square-root T c (p) curve for comparison. T c (p) data from 
Ref . [13] (solid circles) are sparse in this range. Vertical bars on 
the T c data show transition widths estimated from rounding 
of cti (T) fluctuation peaks. The corresponding doping spreads 
(horizontal bars) show little variation with doping. 



tuations near optimal doping [2H US] together sug- 
gest that close to the critical doping, p c , the transi- 
tion out of the superconducting state may be controlled 
by fluctuations near a quantum critical point. In the 
scaling theory of a quantum phase transition, physical 
properties are related to a single, divergent correlation 
length, £ cx |a;| _l/ , where x = (p — p c ). Scaling anal- 
ysis of the XY model predicts T c cx £ -z cx \x\ vz and 
Pso £~( d ~ 2+2: ) cx \x\ u ( d ~ 2+z \ where z is the dynam- 
ical critical exponent [5J ISO]- Together, these rela- 
tions imply T c cx p*Q 2 + z ) _ independent of v. In d = 2 
this requires T c cx p s o, whereas our experiments show 
T c ~ pIq 2 ■ In d = 3, T c cx p% < ' 1+z ' > and is consistent 
with our observations if z = 1, the so-called (3 + 1)D- 
XY universality class [5]. D — 4 is the upper critical 
dimension of the XY model, so mean-field critical be- 
haviour (y — |) should follow. However, this choice 
of v is not compatible with the doping dependences of 
T c and p s0 shown in Fig. |3j scaling arguments predict 
T c cx a; 1 / 2 and p s o cx x whereas the analysis shown in 
Fig. |] has T c cx x and therefore p S Q cx x . This diffi- 
culty may well stem from our determination of doping 
in Fig. [3] which is based on a linear fit and extrapola- 
tion of the T c (p) data from Ref. [131 The sparseness of 
the T c (p) data, along with the lack of an anchor point 
as T c — > 0, allow a different interpretation in the low 
doping regime. In Fig. [3j all data with T c > 8 K fol- 
low an accurately linear doping dependence of p s o, which 
seems to extrapolate well to much higher dopings. The 
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FIG. 5: (color online). Fluctuation peaks in <Ji(T) at 2.64 GHz as doping is varied. Peak width is approximately constant in 
the higher doping range and then broadens considerably below T c w 8 K, consistent with a steepening of T c (p) in that range. 



robustness of this form leads us to put forward the fol- 
lowing suggestion as a means of understanding the data. 
We propose that p s o(p) is, in actual fact, proportional 
to (p — p c ) in the low doping range and use this linear 
mapping to determine T c (j>) . The results of this analysis 
are shown in Fig. [4] The mapping leaves T c (p) consistent 
with the data of Ref. [T3J the effect is simply to refine the 
T c (p) curve in the vicinity of the critical doping. T c (p) 
now grows initially as (p — Pc) 1 / 2 , in accord with scal- 
ing arguments, before crossing over to a linear doping 
dependence at higher dopings. Over a substantial dop- 
ing range this yields the well-known result that T c scales 
approximately with p s o [9 . In addition, our proposed 
interpretation gives a plausible explanation of another 
aspect of the data. At the lowest dopings, the fluctua- 
tion peaks in <J\(T) broaden considerably, as shown in 
Figs. [4] and [5j In the presence of small, macroscopic vari- 
ations in oxygen concentration, such broadening would 
be a natural consequence of a steepening of T c (p) on the 
approach to p c . Future experiments, including those of 
the sort carried out in Ref. [13j will be needed to con- 
firm our conjecture. However, the proposed scenario has 
several compelling features, not least of which is a sub- 
stantial simplification of our theoretical picture of the 
transition from nonsuperconductor to superconductor in 
the underdoped cuprates. 

In summary, at the critical doping for superconduc- 
tivity p S Q becomes nonzero and appears to grow linearly 
with doping, at a rate that remains constant up to much 
higher doping. Fluctuations in the low doped range are 
three dimensional, with no indication of the physics of 
2D vortex unbinding, and with critical exponents char- 
acteristic of (3 + 1)D-XY universality. Quantum fluctu- 
ations appear to control T c in the critical region, with 
T c oc (p — pc) 1 / 2 , but become less effective at depleting 
p s away from p c , allowing T c (p) to cross over to a linear 
doping dependence. 
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